Influence of polymer excluded volume on the phase behavior of colloid— polymer 

mixtures 
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We determine the depletion-induced phase-behavior of hard sphere colloids and interacting poly- 
mers by large-scale Monte Carlo simulations using very accurate coarse-graining techniques. A 
comparison with standard Asakura-Oosawa model theories and simulations shows that including ex- 
cluded volume interactions between polymers leads to qualitative differences in the phase diagrams. 
These effects become increasingly important for larger relative polymer size. Our simulations results 
agree quantitatively with recent experiments. 
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Adding a sufficient amount of non-adsorbing poly- 
mer chains to a stable colloidal dispersion can 
cause a depletion-induced separation of the disper- 
sion into colloid-rich and colloid-poor phases, a strik- 
ing effect which has been extensively investigated 
experimentally ||, 0, |], [|, [| J?], || and which has 
important industrial]^ and biological [fiot implications. 
The entropy-driven depletion attraction between col- 
loidal particles can be tuned by varying the polymer- 
to-colloid size ratio q = R g /R c (where R g is the poly- 
mer radius of gyration and R c the radius of the spherical 
colloids) and the polymer concentration, thus providing 
a unique opportunity of generating rich phase behavior. 
A theoretical description of colloid-polymer mixtures is 
a challenging Statistical Mechanical problem because of 
the large length scale differences between the size of the 
colloids and the polymer segments. The simplest and 
most widely made assumption is to consider ideal (non- 
interacting) polymers obeying Gaussian or random walk 
statistics fill. A further simplification is provided by the 
Asakura-Oosawa (AO) model jl^] whereby polymer coils 
arc treated as mutually penetrable spheres, which are 
excluded from a sphere of radius (R c + R g ) around each 
colloid. Gast et al.|l3|j and Lckkerkerker et al.|l4| used 
this model to calculate the phase diagrams of colloid- 
polymer systems. Computer simulations of hard sphere 
(HS) colloids and ideal lattice polymers p5[ yielded re- 
sults in good agreement with predictions based on the 
simpler AO model. Thus the phase behavior of mixtures 
of HS colloids and ideal polymers is well understood, at 
least for size ratios q < 1. However, polymers rarely be- 
have as ideal, except perhaps near the theta point [pd|. 
The more general problem involving interacting polymer 
chains is much more difficult and over half a century of 
theoretical work has shown that excluded volume inter- 
actions between monomers lead to important qualitative 
and quantitative differences in the properties of polymer 
solutions This letter presents the first large scale 
systematic simulations of the phase behavior of mixtures 
of colloids and interacting polymers for three size ratios 
q, and compares the calculated phase diagrams to recent 
experimental dataj| |], @, ||. 

Some earlier attempts to account for polymer interac- 



tions in colloid-polymer mixtures were based on a pertur- 
bation theory around 0-point conditions fl(| or on integral 
equations p7|. In the present work the self- avoiding walk 
(SAW) model is adopted for the interacting polymers, 
which is known to be a very good representation of poly- 
mers in good solvent ]TH . A full scale simulation of N c 
colloidal hard sphere particles and N p polymer chains, 
each made up of L > 1 monomers or Kuhn segments 
would be a daunting task. However, large-scale simula- 
tions become feasible within a coarse-grained description 
of the polymers, whereby the latter are represented as 
single particles interacting via an effective pair potential 
between their centers of mass (CM). Such effective pair 
potentials can be calculated by a well controlled tracing 
out of the individual monomer degrees of freedom along 
the "polymer as soft colloids" approach we have recently 
put forward g§ |9[ §(j, ||. This approach was shown 
to reproduce, within statistical simulation uncertainties, 
the correct equation of state of pure interacting polymer 
solutions jl^, pH , as well as the correct one-body free en- 
ergy of inserting a single colloid into a polymer solution, 
and the related polymer surface tension p2|. Compared 
to very time-consuming monomer-level simulations, the 
coarse-graining method moreover yields accurate deple- 
tion potentials between two plates or between two 
colloidal spheres |2^| for polymer volume fractions up to 
<f>p = ^xRgNp/V w 2. The success with the two-colloid 
problem suggests that the coarse-graining procedure may 
be fruitfully extended to the full many-body problem of 
the phase behavior of polymer-colloid mixtures as long 
as 4> p is not much greater than 1 (which marks the cross- 
over to the semi-dilute regime). 

Firstly, as described in our previous publications, the 
effective potentials for the polymer-polymer interaction 
were obtained from simulations of a bulk system of self 
avoiding walks (SAW), at various concentrations. All 
simulations were for polymer chains of length L = 500 
segments, with the zero-concentration R g — 16.83 lattice 
units. This length is sufficient to show proper scaling be- 
havior in the semidilute regime pj| . The concentration- 
dependent effective interactions Vpp(r;<f> p ) were obtained 
by inversion of the CM radial distribution function g(r) 
using the hypernetted chain (HNC) integral equation J19[, 
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FIG. 1: Phase diagrams in the polymer reservoir concentration-colloid volume fraction representation. The open symbols 
denote the fluid-fluid and the dashed lines with error bars the fluid-solid coexistence curves for interacting polymers. The 
filled symbols denote the fluid-fluid phase-lines for the two-component AO model. The critical points are indicated by a larger 
open and closed symbol for the interacting polymers and AO model, respectively. The thin solid curves represent free-volume 
perturbation theory phase diagram, with a cross denoting the location of the critical point. The dashed line in the left panel 
is the free volume theoretical metastable fluid-fluid binodal. 



and were accurately parameterized by sums of Gaus- 
sians 

The concentration-dependent potentials v cp {r; <f> p ) for 
the colloid-polymer interaction at each q and <p p were 
obtained from simulations of a single hard sphere in a 
solution of SAW polymers. The CM concentration profile 
was inverted by using a two-component version of the 
HNC equation, and the resulting v C p(r;(f> p ) was fitted to 
an exponential form [ pl| . 

The direct colloid-colloid interaction v cc (r) was taken 
to be hard-sphere like. These three interparticle poten- 
tials are the basis of the coarse-graining scheme. Each 
polymer is now reduced to a single effective particle, 
opening the way to large-scale simulations of a binary 
mixture of polymers and colloids. It is important to note 
that in a system with a finite density of colloids, the poly- 
mer concentration parameter 4> p in the effective poten- 
tials v pp (r;4>p) and v cp {r;<j) p ) must be chosen to be that 
in a reservoir of a pure polymer system in osmotic equilib- 
rium with the two-component system of interest. In other 
words, as discussed in detail by other authors 0, ^4|, p5[ , 
the effective interactions should be taken at the chemical 
potential [i p of the polymers in an osmotic reservoir. 

The Gibbs-Ensemble Monte Carlo (GEMC) 
technique |2?j is naturally suited for studying 
fluid-fluid phase separation. The chemical potential ji v 
was fixed by a standard grand canonical prescription |2q| 
and the number of colloidal particles were fixed at 
N c = 108, N c = 150 and N c = 200 for polymer-colloid 
size ratios of q = 0.34, q = 0.67 and q = 1.05 respec- 
tively. The GEMC simulations yielded histograms for 
the polymer and colloid densities in both boxes. At 
chemical potentials above a critical value, the two boxes 
show different colloid densities, corresponding to phase 



separation. 

Besides the fluid-fluid phase separation, colloid crys- 
tallization can also occur in colloid-polymer mixtures. 
At zero polymer concentration the fluid-solid coexistence 
occurs at colloid volume fractions of <j> c = 0.494 and 
4> c = 0.545, as expected for a pure HS system. The 
effect of the added polymer is initially to widen the co- 
existence curve. At very high polymer concentration, a 
dense colloidal crystal can be in equilibrium with a very 
dilute colloidal "gas" (see e.g. Ref . p9| ) . Experiments 
and previous approximate theoretical work show that if 
q < 0.3 the fluid-fluid critical point becomes metastable 
w.r.t. the crystallization phase-line. For this reason, it 
is important to also calculate the fluid-solid phase lines, 
which was done using Kofke's Gibbs-Duhem integration 
technique |2^, |30|. Starting with zero polymer activity 
at HS fluid-solid coexistence, we performed a series of 
Nfi p PT ensemble simulations, integrating the Clausius- 
Clapeyron equation dP/dfi p — AN p /AV along the co- 
existence line, where AN p and AT^ are the differences 
in number of polymers and in volume between the two 
phases respectively. 

To compare with ideal polymer theories, full GEMC 
simulations were also performed for the AO penetrable 
sphere model of polymers with the same numbers of 
colloids and size-ratios as for the interacting polymer sys- 
tems. These simulations should provide an accurate rep- 
resentation of a true ideal polymer system [T^|. For com- 
parison, we also calculated the phase-diagrams within the 
free- volume theory of Lekkerkerker et al. H] . 

The complete phase diagrams for the ideal and inter- 
acting polymer models at all three size ratio's are plot- 
ted in Fig. [l], in the polymer reservoir concentration- 
colloid volume fraction representation. Although the 
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polymer chemical potential is the natural control variable 
in the simulations, the reservoir polymer concentration 

(r) 

<fip , directly obtainable from \i v through the SAW equa- 
tion of state(ll], is a more suitable variable because 
this brings the interacting and non-interacting polymer 
phase-lines much closer to each other than a direct com- 
parison of the chemical potentials would. 

For q = 0.34 the fluid-fluid binodal has just become 
metastable. This result is consistent with experiments [||, 
and close to the prediction of the AO model and other 
non-interacting polymer theories |l4|, [l^]. The direct AO 
model simulations show a fairly good agreement with the 
interacting polymer simulations and are in good qualita- 
tive agreement with the free- volume perturbation theory, 
although the latter overestimates the critical colloid den- 
sity. Although we did not explicitly calculate the triple 
point for the two-component AO system, it is expected to 
be fairly well located by the simpler free volume theory 
shown in the figure |l4], |l5| |. For q = 0.67 and q = 1.05 
larger differences are seen between the interacting poly- 
mers and the non-interacting AO model system. The 
triple point predicted by the latter moves to much higher 
polymer reservoir concentration than that found for in- 
teracting polymers. Whereas the critical point for inter- 
acting polymers stays near <p c = 0.2, the ideal polymer 
critical point moves to lower colloid densities, an effect 
that was predicted in the literature Jjl]. Finally we point 
out that since the polymer reservoir concentration stays 
well below the concentration <j) p — 2, i.e. in the regime 
where we previously found good agreement for the two- 
body problem, we expect our coarse-graining model to 
provide a very accurate representation of the fully inter- 
acting polymer-colloid system. Should the critical poly- 
mer concentration continue to rise with increasing q, then 
the coarse-graining method would become less trustwor- 
thy for large size ratios. 

Experiments are usually done at fixed polymer con- 
centration, so that the tie-lines are no longer horizontal. 
This representation is shown in Fig for q — 0.67. Triple 
points turn into triangular areas in which three phases co- 
exist. These have been measured experimentally using a 
mixture of PMMA-latex and polystyrene |jj . We compare 
the experimental triple point results for q = 0.57 to the 
theoretical q — 0.67 diagram. Although the systems are 
not entirely equivalent, a much better agreement is found 
with the interacting polymer simulations than with the 
AO model estimates. The simulations also agree quanti- 
tatively with the fluid-fluid binodal which was accurately 
measured for q = 0.677 in recent experiments on silica 
particles in toluene [||. Similar quantitative agreement 
(not shown) was also found between the q — 0.34 simu- 
lations and q = 0.377 experimental binodals||. 

It is also instructive to compare the fluid-fluid bin- 
odals for different size-ratios, as done in Fig. ||. For 
small size ratio the AO and the interacting polymer es- 
timates are close, as also seen in Fig. [l| This is not 
surprising, since we have shown previously that a can- 
cellation of errors in the AO model leads to fairly good 
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FIG. 2: Phase diagram in the polymer volume fraction- 
colloid volume fraction representation for q = 0.67. Symbols 
as in Fig 1. The large triangular areas denote the estimates 
for the triple point: the shaded area corresponds to the exper- 
imental data [7]. The simulated binodal agrees quantitatively 
with the experimental data of ref. h 
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FIG. 3: Fluid-fluid binodals in the polymer volume fraction- 
colloid volume fraction plane. Symbols as in Fig 1. The 
crosses indicate the estimated position of the critical point. 
Note the qualitative differences in the effect of increasing q 
on the critical points as indicated by the arrows. 



effective pair potentials for small q and <ft p p3fl. But the 
simulations for q > 0.34 exhibit qualitatively different 
behavior: the critical colloid density of the AO model 
decreases whereas the polymer concentration does not 
change much. Including polymer interactions has the op- 
posite effect: the critical polymer concentration increases 
and the colloid density hardly changes. Recent integral 
equation calculations^ [ijj of the spinodal phase lines 
also show an increase in the polymer concentration with 
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FIG. 4: The behavior of the critical concentrations as a func- 
tion of size ratio q. The interacting polymer predictions for 
the colloid critical volume fraction 4>i cp ^ agree with experi- 
mentally determined critical concentrations. The free volume 
theory overestimates 4>i cp ^ for small size ratio but underes- 
timates it for large size ratio. The critical polymer volume 



fractions 



for the interacting polymer simulations are in 



reasonable agreement with the experiments of Ref. ^] 
compare less favorably to those of Refs fet] and M] . 



but 



increasing q. 

The behavior of the critical point is summarized in 
Fig. ||, where the predicted trend for the critical col- 
loid density is shown to agree very well with another 
set of experiments^, ||. We have recently performed di- 
rect simulations for much larger size ratios (q > 5), and 
found that the colloid critical point is near the same den- 
sity as found here[^2). This suggests that even in the so 
called "protein limit", q — > oo, the colloid critical den- 
sity is finite. The situation for the critical polymer vol- 



ume fraction is less clear. The present simulations seem 
to agree fairly well with the experimental estimates from 
Ref. [|| fjj although we should stress that the critical point 
was not determined accurately. Two other experimental 
results from Ref. Q| and Ref. || based on mixtures of sil- 
ica and PDMS polymer are also included. These studies 
determined the critical point accurately, and the colloid 
density agrees well with our simulations. However, the 
polymer concentrations are twice as large, which might 
be due to the highly polydisperse polymers used by these 
authors (M w /M n = 2, versus M w /M n = 1.04 in Ref 0). 

In summary, this letter presents the first large scale 
computer simulations of the full equilibrium phase- 
diagram of a binary mixture of colloids and interacting 
polymers. These simulations would not have been pos- 
sible without using the accurate "polymers as soft col- 
loids" coarse-graining approach. Even though the phase- 
separation happens mainly in the dilute regime of the 
polymer solution, we find important qualitative differ- 
ences with ideal polymer behavior. In particular, the ab- 
solute polymer concentration <f> p at the critical point in- 
creases with increasing q for interacting polymers, while 
it decreases with increasing q for non-interacting poly- 
mers. Similarly, the critical colloid packing fraction 
remains nearly constant for interacting polymers, while it 
decreases with increasing q for non- interacting polymers. 
We also showed that by including excluded volume inter- 
actions wc find quantitative agreement with experiments. 
In conclusion then, just as excluded volume interactions 
are known to strongly influence the behavior of polymer 
solutions, we find here that polymer excluded volume 
interactions have important qualitative effects on the be- 
havior of polymer-colloid mixtures. 
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